1 



Advances in the Determination of Quark Masses 

T. Bhattacharya and R. Gupta a * 

a MS B285, T-8, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA. 

1 Significant progress has been made in the determination of the light quark masses, using both lattice QCD and 

sum rule methods, in the last year. We discuss the different methods and review the status of current results. 
0^ ■ Finally, we review the calculation of bottom and charm masses. 



o 

Q 

m 

(N 
> 
in 
a\ 
o 
o 

On 



X 



1. Introduction 

One of the primary goals of phenomenology 
is to determine the fundamental parameters of 
the standard model. Of these, the five quark 
masses m u , mj, m s , m c , to& are amongst the least 
well known, and cannot be measured in experi- 
ments. They have to be inferred from the pattern 
of the observed hadron spectrum (chiral pertur- 
bation theory (xPT), HQET, and lattice QCD 
approaches) or from the study of 2-point corre- 
lation functions (sum rules and lattice QCD ap- 
proaches). This review is mainly an assessment of 
lattice QCD results, however we shall also sum- 
marize the results obtained using xPT and sum- 
rules. In particular we shall evaluate the extent 
to which these methods agree with each other and 
discuss recent developments in each. 

%PT relates pseudoscalar meson masses to m u , 
rrid, and m s . However, due to the presence of an 
overall unknown scale in the chiral Lagrangian, 
%PT can predict only two ratios amongst the 
three light quark masses & 



2m s /(m u + rrid) 
m u /md 



Lowest order 
24.2 - 25.9 
0.55 



Next order 
24.4(1.5) 
0.553(43) . 



These ratios have been calculated neglecting 
the Kaplan-Manohar symmetry jj]. This is jus- 
tified by assuming that the higher order terms 
are small as discussed in ||. These ratios, 
when combined with an absolute value deter- 
mined from sum-rules completed the "standard" 

*Editors' note: T. Bhattacharya was originally invited to 
give this talk at LATTICE 97, but could not attend. He 
was replaced by R. Gupta. They agreed to a joint writeup. 



scenario. The status of sum-rule estimates, as 
of 1996 was, 2m = m u + m d = 9.4(1.8)MeV §] 
and m s = 126(13)MeV §, evaluated in the MS 
scheme at scale \i = 2GeV. These, combined with 
xPT ratios, give the standard estimates (in MeV) 



Input m u 
m s 3.7(1.0) 
m 3.3(1.3) 



m d 
6.7(0.8) 
6.1(1.1) 



m 

5.2(0.6) 



115(23) 



These values are consistent with observed spec- 
tra, i.e. one can explain the electromagnetic 
splittings, and the SU(3) breaking in the mesons 
and in the baryon octet and decuplet without 
recourse to large non- linear terms in m q . The 
limitation of this combined analysis, especially 
of the absolute numbers, is that there are no 
independent checks as the quark masses enter in 
phenomenology only in combination with other 
unknown quantities like the quark condensate. In 
the last year there has been considerable activity 
in sum-rule analyses, and we shall summarize 
these in Section 0. 

Lattice QCD is a relative newcomer. Last 
year a major step forward was taken in the un- 
derstanding and quantification of systematic er- 
rors. An analysis of the global data consolidated 
the lattice predictions of significantly lower light 
quark masses 

Quenched Nf — 2 
(m u + m d )/2 3.4(7) ~ 2.7 MeV 
m s 100(31) - 70 MeV. 

The lattice results for the quenched theory were 
considered reliable as three different discretiza- 
tions of the Dirac action, Wilson, clover, and 
staggered, gave results consistent after extrapo- 
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lation to the continuum limit. The unquenched 
estimates (Nt = 2 theory results) were prelimi- 
nary. The main message was that including two 
flavors of dynamical quarks lowers the quenched 
estimates by ~ 20%. 

The abbreviation "1996 data" will be used to 
denote lattice data analyzed in 0. This will be 
used as the benchmark against which progress 
will be measured. Other reviews and references 
to lattice results can be found in . The main 

sources of uncertainty in the lattice estimates we 
shall explore are (i) the extrapolation of the data 
to the continuum limit, (ii) the matching con- 
stants between the lattice and continuum scheme, 
and (iii) the effect of sea quarks (corrections due 
to quenching) . We are happy to report that there 
has been significant progress in all three areas in 
the last year. 

2. Lattice Approach: From Spectroscopy 

Lattice calculations need to determine four in- 
dependent quantities to predict the three light 
quark masses. The fourth is needed to fix the 
lattice spacing a. To distinguish between m u and 
nid, it is necessary to include electromagnetic cor- 
rections. This can be done by including a U(l) 
field in the simulations [Q, however, since the 
effect is small compared to statistical and sys- 
tematic errors, current lattice simulations have, 
for the most part, neglected it. Thus, this review 
is restricted to the determination of the isospin 
symmetric combination to. 

A brief outline of how the three quantities to, 
to s , and a are determined from the light hadron 
spectrum is as follows. Using xPT as the guid- 
ing principle, one writes down the most general 
expansion of hadron masses in terms of quark 
masses toi , 777,2 , 771.3 > 



Mia 2 
M p a 

M s +a 
M A a 



Art + B n a(mi + m 2 )/2 + . . . 
A p + B p a{mx + m 2 )/2 + . . . 
A s + + \Fam x + 2(F - D)am 3 
Aa + BAa(mi + m 2 + to, 3 )/3 + 



,(1) 



where for brevity we shall use 7r, p, A to denote 
members of the pseudoscalar and vector octet and 
the baryon decuplet respectively. Chiral symme- 



try implies A^ = 0. We leave it as a free pa- 
rameter in the fits - it gives a measure of the 
uncertainty in the zero of the mass scale. 

The simplest scenario is that there are no 
higher order corrections to Eq. |l|. Then, 
any triplet like B n ,A p ,B p (or equivalently 
-Btt , A a , -Ba ) can be used to determine to, to s , 
and a. For example, using the first set, the quark 
masses are determined as 
_ Ml _ 2(M K , 



Apa- 1 ) 



Bird 



B„ 



(2) 



and the scale from M p by solving the quadratic 



.1, 



MJ- 



.1 



B n 



Ml 



0. 



(•3) 



A p (-) 2 - -., 
a a H-n 

Similar relations exist for other choices of observ- 
ables. The scale 1/a can, in fact, be taken from 
other observables like M n , Ma, fn, string ten- 
sion, ro, or the IP— IS* splitting in quarkonia. We 
choose a(M p ) based on its ready availability and 
statistical quality of M p versus M n , Ma and /„• 
data. Different choices lead to different results, an 
unavoidable uncertainty inherent in the quenched 
approximation. We shall refer to this hadron 
spectroscopy method by the abbreviation HS. 

If linearity were exact, or the exact expansions 
in Eq. [ljwere known, then any single set, like the 
pseudoscalar octet masses, could be used to fix 
both ma and m s a. There is no a priori reason 
to assume that the higher order chiral corrections 
are negligible. In fact, using M^/Mk or M v /Mk 
give different estimates for m s /m in lowest or- 
der xPT. Also, present lattice data show non- 
linearities in the pseudoscalar and vector meson 
data if a sufficiently large range of quark masses is 
chosen [ ^3pl| ]. However, these non-linearities are 
small, and are very hard to detect in the typical 
range m s /3 < m q < 2m s used to extract quark 
masses. Thus, unless otherwise stated, we have 
used just linear fits in the chiral extrapolation. 

The more serious problem facing quenched sim- 
ulations is that this range cannot be extended 
easily to much smaller quark masses due to the 
presence of artifacts called quenched chiral logs. 
(For recent reviews and references to this body 
of work see 0,^5).) These terms, which are not 
present in the normal chiral expansion, are singu- 
lar in the limit m q — > and are a consequence of 
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the fact that, in the quenched approximation, the 
rj' propagator has a single and a double pole. The 
approach, therefore, has been to fit the quenched 
data in this limited range, where these artifacts 
are small, keeping just the normal chiral expan- 
sion. The quenching error is then the change in 
these coefficients as sea quark effects are added. 

Finally, a test of the combined uncertainties 
due to quenching, and chiral and a — > extrapo- 
lations is that quark masses extracted using dif- 
ferent observables should all give the same re- 
sults. We shall discuss the extent to which this 
is satisfied by comparing m s extracted from Mk 
with that from Mk* or (labeled m s (Mjf), 
m,(Mx*), m s (M,p) respectively). Most of the fig- 
ures we shall present will be for to, however, note 
that due to the linear approximation in Eq. [j], the 
result for m s (Mj() is = 25.9m. m s (M^) gives an 
independent estimate, but the statistical signal in 
the data for vector state masses is not as good as 
for pseudoscalars. This will be evident from the 
figures we show for m s (M^). 

2.1. Definition of quark mass 

For staggered fermions, the fits in Eq. [I] are 
made as a function of the bare quark mass input 
into the simulations. For Wilson-like fermions, 
we define light quark masses by 



am bare = log(l + (—- -!-)) 



(4) 



where k c is the critical value of the hopping pa- 
rameter at which the lattice pion mass vanishes. 
(Reference [Q used ambare = 1/2/c — 1/2k c , which 
is consistent to 0(a). This change leads to differ- 
ences of a few percent.) So, compared to stag- 
gered fermions, Wilson like fermions require the 
determination of an additional quantity n c . The 
need to calculate k c can be avoided by using the 
Ward Identity method described below. From 
mbare, the MS mass is obtained as 



/JI- 



MS 



Z m m ba 



(5) 



where Z m = 1/Zg is the matching factor between 
MS and lattice schemes. Since this factor turns 
out to be large at 1-loop, its reliability has to be 
checked by non-perturbative methods. 



2.2. Ward Identity (WI) method 

For Wilson like fermions, one can extract quark 
masses from pseudoscalar mesons by using the 
axial vector ward identity 



ZaO^A^x) = (mi + m 2 )Z P P(x) + 



(6) 



where is the appropriate flavor non-degenerate 
bare axial current ip^p^fp, P is the corresponding 
pseudoscalar density ip^ipy an< 4 the Z's are the 
corresponding renormalization constants. The 
dots represent discretization corrections, whose 
size depends on the order to which improvement 
of the action and operators has been carried out. 
The renormalized quark masses are then given by 
the ratio of 2-point functions 



(mi 



m 2 ) 



Z A (dMx)J(O)) 
Z P (P(x)J(0)) 



(7) 



where J is any operator that couples to pions. 
The advantage of this WI method is that k c does 
not enter into the calculation. The limitation is 
that only the pseudoscalar sector is tested. 

3. Details of the analysis of World Data 

Different groups have their own favorite ways of 
doing the fits, and of dealing with the two major 
sources of systematic errors, the extrapolation to 
the continuum limit and the determination of Z's. 
Consequently, the same data analyzed by differ- 
ent groups can lead to slightly different estimates 
of quark masses. Since we wish to do a global 
analysis, we have attempted to minimize these 
relative differences. We start with the data for 
pseudoscalar and vector mesons as a function of 
quark masses and g 2 , and redo the analysis to ex- 
tract Ayr, BnjAp, Bp. For the physical masses we 
use M ff = 135 MeV, M p = 770 MeV, M K = 495 
MeV, M K * = 894 MeV, and M = 1020 MeV. 
(Ref. 0| used = 137, thus to quoted there is 
about 3% higher.) This still leaves two sources of 
uncertainty, (i) The difference in lattice sizes and 
the subjective bias in fits made to extract hadron 
masses. This we check by comparing data from 
different collaborations, (ii) The statistical cor- 
relations in the data. These require knowledge 
of the covariance matrix which is not available in 
most cases. 
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To convert the lattice mass to a continuum 
scheme like MS, we require the calculation of 
either Z$ in both lattice and continuum renor- 
malization scheme (HS method), or of Za and 
Zp (WI method). For these matching factors 1- 
loop results have been used in the past. Now, 
non-perturbative estimates are becoming com- 
mon. The calculation of the Lepage-Mackenzie 
a s , "horizontal" matching between the lattice 
and continuum schemes in the perturbative ap- 
proach, details of the error analysis, and the 
evolution in the continuum are the same as in 
. The relevant formula for the 2-loop evolution 
needed with the 1-loop matching is jl6j 

m{Q) fg 2 (Q)Y° /m ( 1] 9 2 (Q)-f(n) A (g) 



16tt 2 



where J = (7i/3o — Yo/3i)/2/3q. From this the 
renormalization group invariant mass fh is de- 
fined as 



fh = m([x) 



2/3o3 2 (^)r 7o/2/? ° f. 3 2 (M) 



16tt 2 



1 



16tt 2 



J .(9) 



This evolution in the continuum or the conversion 
to fh introduce no new lattice uncertainty. 

Now we briefly highlight two important devel- 
opments by the APE and ALPHA collaborations 
that overcome the uncertainty introduced by the 
perturbative matching. 

The APE collaboration has carried through 
non-perturbative determination of the lattice Z's 

fusing the chiral Ward identity method |lS| , 
At present non-perturbative estimates of all 
the necessary lattice Z's for all the actions we dis- 
cuss are not known. As these become available, 
the uncertainty in using the 1-loop relations will 
be removed. We shall discuss the impact of the 
known Z's on current data. 

The ALPHA Collaboration has initiated a very 
successfully non-perturbative program to improve 
the discretization of the Dirac action, fermionic 
operators, and their renormalization constants 
pcfl . The improvement they propose over APE's 
chiral Ward Identity method for determining the 
Z's is to eliminate making any connection with 
a continuum scheme by directly computing the 
renormalization group independent quark mass. 



The essential steps in their method are (i) calcu- 
late the renormalized quark mass using the WI 
method at some infrared lattice scale (Zp and 
Za are also calculated non-perturbatively in the 
Schrodinger functional scheme), (ii) evolve this 
result to some very high scale using the step scal- 
ing function calculated non-perturbatively, and 
(iii) at this high scale, where a s is small, make 
contact with lattice perturbation theory to de- 
fine fh. This non-perturbatively calculated fh is 
scheme independent, thus subsequent conversion 
to a running mass in some scheme and at some 
desired scale can be done in the continuum us- 
ing the most accurate versions of the scale evolu- 
tion equations pl| . No data for fh has yet been 
released. The status of this approach has been 
reviewed by Luscher at this conference pcfl . 

4. Quenched Wilson fermion results 

The Wilson formulation of the Dirac action is 
the simplest and has been the most commonly 
used in numerical simulations until very recently. 
Its disadvantage is that discretization errors and 
violations of chiral symmetry begin at 0(a). Nev- 
ertheless, a 1996 compilation of the world data for 
quark masses |7| showed that an extrapolation to 
the continuum limit can be made keeping only 
the lowest order corrections, giving 



m 



3.4(7)MeV[l - 
94(27)MeV[l 



1.3(2)GeV a] 
- 1.9(7)GeV a}. (10) 



The surprise was the size of the 0(a) errors. The 
typical size of the slope in a obtained in the ex- 
trapolation of other observables like M pi f^, . . . is 
a few hundred MeV, so the 1.2 — 1.7GeV values 
needed to be understood. 

Significant progress in this sector has been 
made by the recent high statistics, large lattice 
calculations by CP-PACS which have been re- 
viewed by T. Yoshie in ]22^] . They have provided 
four new data points on lattices of size > 3 fermi. 
The new data are shown in Fig. |l|. To high- 
light the statistical improvement we show data 
at (3 = 6.0 from the next best calculation (with 
respect to both statistics and lattice size) |l| . 
The data from both HS and WI methods are well 
fit by just a linear correction, and extrapolate to 
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0.2 0.4 0.6 

a (GeV -1 ) 



Figure 1. Linear extrapolation of to versus a(M p ) 
for Wilson fermions using HS and WI methods. 
The 1996 data and fit are shown in green. 

roughly the same value 

to = 4.1(l)MeV[l + 0.49GeV a] HS 

to = 4.2(l)MeV[l - 0.32GeV a] WI. (11) 

The CP-PACS data, therefore, gives a signifi- 
cantly smaller slope and a larger value for to. 

In Fig. [l] we also show the 1996 fit. The change 
in the slope is due to a pivot about j3 « 5.9. Fur- 
ther analysis based on the CP-PACS data shows 
that the three APE points at j3 = 6.3, 6.4, that 
biased the 1996 fits towards lower values of to, 
are not compatible within errors with CP-PACS; 
the estimates of both B v and a differ. Our con- 
clusion is that these data should not be used in 
"modern analyses" . 

The CP-PACS collaboration also calculate m s 
using both Mk and (M^ and Mk* give con- 
sistent results). The new data are shown in Fig. |[ 
The data show that to s (tox) and m s {M§) are dif- 
ferent even after extrapolation to the continuum 
limit, as was the case in the 1996 data, 

m s (M<f,) = 139(ll)MeV[l + 0.47GeV a] 
m s (M K ) = 107(2)MeV[l + 0.49GeV a}. (12) 

Assuming that this difference is not an artifact of 
statistical errors or due to the uncertainty in the 
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Figure 2. Linear extrapolation of m s (M^) and 
m s (M ( f > ) versus a(M p ) for Wilson fermions. The 
1996 fits are shown in green for comparison. 

extrapolation, it could be due either to quench- 
ing, or to the failure of the assumption of linearity 
in the chiral fits. In that case one will need un- 
quenched data to resolve this issue. Meanwhile, 
we consider this difference as one of the remain- 
ing systematic errors. These estimates for m s are 
larger than those reported in 1996. The reasons 
for the change are the same as for to. 

The second important contribution of the CP- 
PACS analysis is the reconciliation of the HS 
and WI methods. The data, and the fits shown 
in Fig. |l| indicate that the difference between 
the two methods at finite a is due to different 
O(a) discretization errors. Giusti et al. JlTj ar- 
gue that this difference is due to the failure of 
the 1-loop Z's. We have therefore reanalyzed 
the CP-PACS/LANL WI data using the non- 
perturbative Za,Zp calculated by Giusti et al. 
(we interpolated or extrapolated the Giusti et al. 
results to other f3 linearly). These reweighted 
points are plotted using the symbol octagon in 
Fig. [|. What is remarkable is the agreement be- 
tween non-perturbative WI and perturbative HS 
data around /3 = 6.1, and the change in the slope, 
i.e. the slope switches sign and ends up being even 
larger than that for HS with 1-loop Z's! This 
equality of HS and WI data around (3 = 6.1 is con- 
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sistent with the findi ngs of Giusti et al. ||17(] as dis- 
cussed in section 4.3. What the APE data cannot 



expose is the slope in a from simulations at just 
ft = 6.0 and 6.2. We would also like to point out 
that z no "-P ert /ZP ert typically decreases rapidly 
with ft as, for example, in Table Hence, if 
the 1-loop Z m turns out to be an underestimate, 
then correcting for it will increase the slope of 
the HS data as well. Thus, the issue of the size 
of the slope, and of the final extrapolated value 
will be resolved soon, once the non-perturbative 
estimates of Z m are made available pij ]. 

4.1. Clover Action 

One way to reduce/remove the large 0(a) cor- 
rections in the Wilson formulation, and thus im- 
prove the reliability of the a = extrapola- 
tion, is to simulate the Sheikholeslami-Wohlert 
(SW or clover) action. In the last couple of 
years a number of calculations have been done 
using this action. Unfortunately, different calcu- 
lations use different values of Csw- (i) tree-level 
value Csw = 1, (h) tree- level Tadpole Improved 
(TI), (iii) 1-loop tadpole improved, and (iv) non- 
perturbative O(a) improved. To get a feel for the 
effects of tuning Csw > we first show in Fig. || the 
data for m for different values of Csw as a func- 
tion of a. The data at ft = 6.0 show that as Csw 
is increased, m decreases and a, as determined 
from M p , increases. The expectation is that as 
(3 — > oo, the ordering should stay the same, only 
the spread should decrease. 

The most extensive data are using the tadpole 
improved clover action (we do not distinguish 
between tree-level and 1-loop improved as the 
difference in Csw is small and we assume that it 
amounts to a negligible change in quark masses). 
These data for m from JLQCD ||,||, LANL 
f25|| , and UKQCD |2(| collaborations are shown 
in Fig. [4. We also show the Fermilab data 
given in Of the « 10% difference between 
JLQCD and Fermilab data at (3 = 5.9 and 6.1, 
k, 6% comes from the setting of a (Fermilab uses 
IP — IS splitting in charmonium) , and the rest 
from differences in B n and conversion to MS. 
Since the data needed for extracting a(M p ) for 
the Fermilab runs are not available, we cannot 
use these points in our continuum extrapolation. 
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Figure 3. Behavior of m and a(M p ) as a function 
of the clover coefficient Csw at ft — 6.0. Also 
shown are the CP-PACS data and fit from Fig. |. 

A fit to the rest of the data, assuming that the 
errors are 0(aa), gives 



m = 3.84(10)MeV[l + l.OGeV a(a) a] 



(13) 



shown by the fancy square at a — 0. An analysis 
of m s (MK) and m s {M ( j > ) gives 

m s (M ) = 117(8)MeV[l + 1.5GeVa(a)a] 
m s {M K ) = 99(3)MeV[l + 1.0GeVa(o)o]. (14) 

In the figure we plot, for simplicity, the linear fit 
which has roughly the same x 2 



■ = 3.72(13)MeV[l + 0.26GeV 



(15) 



This lies ~ lc below our preferred fit and value 
given in Eq. [13]. Similar linear fits to the m s (MK) 
and m s (M c f,) data are shown in Fig. [|. 

The APETO ||, QCDSF @ and UKQCD 
|p6|| collaborations have calculated quark masses 
using the non-perturbative value of Csw deter- 
mined by the ALPHA collaboration. The data 
for m by QCDSF and UKQCD at ft = 6.0 and 
6.2 are shown in Fig. ^ and compared to tadpole 
improved clover. We find the following features. 
At ft = 6.0, data with the WI (QCDSF) and HS 
(larger volume QCDSF and UKQCD) methods 
agree, and roughly lie on the tadpole improved 
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Figure 4. m from the TI clover data versus 
a(M p ). Also shown are CP-PACS Wilson data 
and fit from Fig. [l]. Result of extrapolation in aa 
is shown by the fancy square at a = 0. 



clover line. Going to (3 = 6.2, the QCDSF WI 
data show a small increase. Consequently, a fit 
to the two QCDSF- WI points, assuming 0(a 2 ) 
errors, gives the large value m = 5.1(2) MeV (the 
value in the Fig. ^| is slightly different as it is from 
our analysis). On the other hand the QCDSF-HS 
value decreases slightly, while the UKQCD-HS 
data (still preliminary) show no a dependence. 
Our overall conclusion is that the range in j3 is 
too narrow to allow a meaningful extrapolation 
in a with just two points. 

Both APETO |7| and QCDSF |§ have ex- 
tracted tti s at (3 — 6.2 on 24 3 x 48 lattices. Their 
raw data for pseudoscalar and vector masses are 
consistent. Since the extraction of m s includes 
ma corrections to the renormalization constants, 
their results are correct to O(a). The estimates 
of m s are 111(15) and 102(2) MeV respectively. 
QCDSF find a small increase between /3 = 6.0 and 
(3 = 6.2, and their value extrapolated to a = is 
112(5) MeV. Higher statistics and larger lattices 
data at other values of f3 are expected in the next 
year. These will improve the reliability of the ex- 
trapolation to the continuum limit. 
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Figure 5. Data for m s (Mic) and m s (M t f > ) versus 
a(Mp) for TI clover fermions, along with simple 
linear extrapolation. Also shown are APE data 
with Csw = 1 for comparison. 

4.2. Staggered Fermions 

The advantage of staggered fermion formula- 
tion is the remnant chiral symmetry that guar- 
antees Ward identities as in the continuum. As 
a result Za = 1 and Zg — Zp, and the quark 
mass is only multiplicatively renormalized. Con- 
sequently, the two methods for determining quark 
masses, HS and WI, are the same. The sticky 
point is that the finite piece in TI 1-loop expres- 
sion Z m = l/Z l ]S cal = 1 - f£(81og0ua) - 39.1) 
is very large. Therefore, the 1-loop matching be- 
tween lattice and continuum may not be reliable. 
Second, the 1996 data at j3 > 6.3 was preliminary 
and suggested a small increase with (3. There 
has been progress on both fronts in the last year. 
First, a partially non-perturbative analysis of the 
reliability of Zp has been done by Gupta, Kilcup, 
and Sharpe pl| ], and secondly there is new data 
by the JLQCD @ and MILC |3l) collaborations. 

The basic idea of the partially non-perturbative 
analysis is that if Z l p cal is large, then one could 
get a more reliable estimate using 



Z 



local 



zsmeared 



2^local 
^smeared 



(16) 



where "smeared" is any discretization of the pseu- 
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Figure 6. to from non-perturbative SW fermions. 
TI clover fit from Fig. is also shown. 
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Figure 7. Staggered fermion data for to versus 
a(m p ). The fit assumes a 2 and a 4 corrections. 



doscalar density that has a reliable perturbative 
expansion. The ratio in the parenthesis, which is 
a large factor, is calculated non-perturbatively. 
On the basis of such an analysis, Gupta, Kilcup, 
and Sharpe found that the 1-loop perturbative 
expression for Z l ° cal is ~ 5% smaller than the 
partially non-perturbative result at (3 = 6.0. 
The caveat in this calculation is that the non- 
perturbative determination of the ratio shows 
large 0(a 2 ) discretization errors, and the extrap- 
olation to a = is based on only two points at 
(3 — 6.0 and 6.2. Thus, the calculation of the 
ratio needs corroboration. For the moment we 
apply this shift when presenting final estimates 
in Section [|. 

The staggered fermion world data are shown in 
Fig. 0. The new data by JLQCD and MILC col- 
laborations are consistent with the old, and con- 
firm the small rise at weaker couplings. Assum- 
ing that the partially non-perturbative analysis 
of Zp validates the 1-loop result, and thus the 
non-monotonic behavior, we fit the data includ- 
ing a 2 and a 4 corrections. (Note that the absence 
of 0(a) corrections implies that the fit approaches 
a = with zero slope.) The results, in MeV, are 

m = 3.35(7) [l-(.7GeVa) 2 +(lGeV a) 4 ] 
m s {M,f,) = 104(5)[l-(lGeVa) 2 + (lGeVa) 4 ](17) 



The discretization corrections are small, and ac- 
count for the non-monotonic behavior, however, 
to resolve this issue a fully non-perturbative cal- 
culation of Z m is needed. The net upshot is that 
the staggered results using 1-loop Z m are basi- 
cally unchanged, and raised by ~ 5% by using 
the partially non-perturbative calculation of Z m . 

4.3. Non-perturbative Z's: APE data 

The APE collaboration has presented new data 
on the estimates of to, to s , and m c using both the 
HS and WI methods p7y . They also provide a 
non-perturbative determination of Za, Zp needed 
in the WI method, however, the Z m used in the 
HS method is still 1-loop. In their calculations 
of the NP Z's for Wilson fermions bare operators 
are used, while for clover fermions a field rotation 
by 1 + ap/4 is included. 

The APE data are summarized in Table [I]. We 
have also included a comparison of the tadpole 
improved Za/Zp with the non-perturbative fac- 
tor (1 + a s C^ n )Z A /Z^ with fi 2 a 2 = 0.9648. 
In keeping with the authors, who consider the 
(3 = 6.4 lattices too small, we neglect the cor- 
responding quark mass data. The first remark- 
able feature is that Wilson and clover fermions 
give consistent results for the WI method. Next, 
based on the consistency of the data at f3 = 6.0 
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Table 1 

APE data for m s (M^) for Wilson (upper box) 
and C'sw = 1 clover fermions. The perturbative 
and non-perturbative Z's are also compared. 





/3 = 6.0 


(3 = 6.2 


(3 = 6.4 


Z m ( pert) 


1.252 


1.222 


1.200 


Z A /Z P { pert) 


1.178 


1.156 


1.140 


Z A /Zp(Npert) 


1.659 


1.507 


1.345 


m s (WI) (MeV) 


128(6) 


117(8) 


107(8) 


m s (HS) (MeV) 


133(16) 


130(16) 


120(16) 


Z m ( pert) 


1.177 


1.157 


1.143 


Za/Z p ( pert) 


1.825 


1.680 


1.611 


Z A /Z P (Npert) 


2.350 


2.033 


1.700 


m s (WI) (MeV) 


125(5) 


127(9) 


106(8) 


m s (HS) (MeV) 


117(7) 


120(9) 


103(9) 



and 6.2, the authors assume that there are no 
significant discretization errors. Thus, their fi- 
nal number m s (MK ) = 123±4±15 is from the 
WI method, averaged over the Wilson and Clover 
fermions data at (3 = 6.0 and 6.2. 

We consider the agreement between WI and 
the HS data fortuitous. If the difference be- 
tween the perturbative and non-perturbative Z m 
is similar to that in Zp, then the final HS re- 
sults would be significantly different. Second, the 
CP-PACS/LANL data show a significant a de- 
pendence which the APE data are not precise nor 
extensive enough to expose. This is the reason 
why their final value is higher than the world av- 
erage given next. 

5. Bottom line on quenched results 

A summary of the quenched results, in MeV, is 



Wilson 

m 4.i(i) 

m s (M K ) 107(2) 
m s (M ) 139(11) 



TI Clover 
3.8(1) 
99(3) 
117(8) 



Staggered 
3.5(1) 
91(2) 
109(5) . 



These differences between Wilson, TI clover, 
and staggered results could easily be accounted 
for by the uncertainties in the extrapolations, 
and/or by the ratio of the non-perturbative to 1- 
loop estimates of Z's. Thus, for our best estimate 
we take the mean and use the spread as the er- 
ror. In case of m s we have an addition uncertainty 
coming from the dependence on the state used to 



extract it, Mr versus Mk> (or equivalently M$). 
This could be due to the quenched approximation 
or an artifact of having used linear chiral extrap- 
olations. Since we do not have control over either 
of these two uncertainties, we again average over 
all the data. To these, we add a second uncer- 
tainty of 10% as due to the determination of the 
scale 1/a. Thus, our best estimates of MS masses 
evaluated at 2 GeV are 



m 
m. 



3.8(4)(4) MeV 
110(20)(11) MeV. 



(18) 



6. N f = 2 Wilson results 



The SESAM collaboration has presented new 
data for n/ = 2 Wilson fermions at (3 = 5.6 
P2j| . The novel feature of their data is that 
pseudoscalar and vector masses have been calcu- 
lated for degenerate (k v = n sea ), and for non- 
degenerate (k v 7^ K sea ) combinations. The su- 
perscripts V and sea refer to valence and sea 
quarks respectively. Their data can be fit to 



Mi 



M n 



A r 



1 



2k v 
1 



Bi 



1 



2K sea 
1 



j^sea 

P 2iK sea 



(19) 



As a result they study the dependence on both 
k v and K sea . 

The first surprise of their calculation is a large 
dependence on K sea (B v ~ B s ). This is sur- 
prising because QCD perturbation theory and the 
success of lowest order xPT suggests that the de- 
pendence on sea quark masses should be small. 

Their second important contribution is that 
they correctly point out that all previous extrac- 
tions of m s from nj = 2 data are incorrect. Pre- 
vious calculations were essentially calculating m s 
in a sea of strange quarks. 

The last issue they raise is specific to Wilson- 
like fermions and concerns the zero of mass scale 
(k c ), and concomitantly the definition of the 
quark mass. They discuss two possible ways of 
analyzing the data, (i) Degenerate extrapolation: 
determine kP c by extrapolating pion masses for the 
degenerate case k v — K sea and measure all quark 
masses with respect to this k°. (ii) Partially 
quenched extrapolation: calculate m, m s and a 
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for a fixed sea quark mass, and then extrapolate 
in the sea quark mass. They argue that the first 
method is the correct one and discard the second. 
It is at this point we disagree with them. What 
we now show is that if one uses the conventional 
perturbative Z m to relate the lattice mass to MS 
scheme, then it is actually the partially quenched 
method that is more appropriate. We also argue 
that if Z m is calculated non-perturbatively, then 
both methods would give the same result for m. 
Lastly, our analysis shows that the large depen- 
dence on K sea is actually that of k c on n sea , and 
not of physical masses. 

The pole mass defined by the inverse propaga- 
tor S^ 1 — Z^(f — m) + • • •, can be written as 

' ' (21)) 



1 o c ( 1 

ma = - — o + dm a 



2k 1 



2k v 2k s 



(21) 



The linear divergence in dm is absorbed in the 
definition of k° (SESAM calls it a ) as 

Then, suppressing the dependence on cx t 

1 1,-/1 1 
ma — — T7 7t + om — 77, 

2k v 2k° c \2n v, 2n sea 

(22) 

Near a = 0, k v and K sea approach i.e. the va- 
lence and the sea quark masses measured in lat- 
tice units approach zero, so we expand 8m about 

,„sea „() 



1 

2^y 



1t& 



rylat 

— A-rr, 



1 

2 K sea 
1 

2K V ~ 2^ 



1 



+ c 



2k v 
1 

2^ 



1 

2^ 



1 



(23) 



where 

Zlat 
m 

1 



i + C 
1 c s 

2^5 



1 



5?)- (24) 



1 + C V \2K sea 

Note that both Z 1 ^ and k c are the partially 
quenched ones. What has happened is that, for 



Wilson like fermions, along with the linear diver- 
gence which is absorbed in 1/2k°, there is a finite 
piece of O(l) that shifts the quark mass. As this 
finite piece is of 0(m) after extrapolation to the 
physical sea quark masses, its effect on the de- 
termination of m is dramatic, while it is a small 
correction in m s . 

In the degenerate case Eq. becomes 



ydeg I J \_ 

' ' 2k 2kI 



(25) 



where Z^ 9 = 1 + ( v + ( sea has a contribution 
from valence and sea quarks. At 1-loop, C sea = 0, 
therefore Z l £ = Z% 9 and k c = k°. What the 
SESAM data are telling us is that corrections to 
this 1-loop result are large, i.e. ( sea ~ 1. Thus, 
we arrive at the conclusion that combining the 1- 
loop Z l £f , which is independent of C seo , with the 
partially quenched extrapolation, which already 
incorporates £ sea as shown in Eq. is more 
appropriate. Another way of saying this is that 
had Z^~ 9 been calculated non-perturbatively, one 
would have found Z^ 9 ~ 2Z^', which would 
have compensated for the smaller m obtained 
from the degenerate extrapolation. 

SESAM, by measuring quark masses with re- 
spect to K% extract m = 2.7(2) MeV and m s = 
140(20). The large shift in m s /m from xPT 
values is attributed to the effect of sea quarks, 
which could change as a — » 0. Our proposal for 
using the partially quenched extrapolation gives 
rfi = 4.7(1) MeV which is ~ 10% below the 
quenched value at the same scale a and roughly 
preserves the xPT ratio. Note, however, when 
using the partially quenched method there ex- 
ist enhanced chiral logs, similar to those in the 
quenched approximation p3| . These afflict the 
theory at small quark masses, but for present 
masses this effect is expected to be negligible. 

If the meson masses are completely indepen- 
dent of K sea , then 



Bi 



/ B V =B ™/ B v, 



V 



(26) 



The SESAM data bears this out within errors. 
However, when this constraint is used in the fits 
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for the vector masses, they find high \ 2 s of about 
50 for 25 degrees of freedom. We take this sta- 
tistically significant, but small, effect to indicate 
that higher order terms in the chiral expansion do 
bring in small dependence of the meson masses on 
the sea quark mass. 

In view of the above discussion, the only n/ = 2 
numbers that survive the above discussed ambi- 
guity from the 1996 analysis are for staggered 
fermions, m ~ 2.7 MeV. More data are necessary 
to establish the continuum limit in that case. For 
Wilson like fermions the final word on the val- 
ues of to and to s has to await data using the WI 
method along with a non-perturbative determi- 
nation of the Z's. 

7. Sum rule determinations of to and m s 

Progress in the sum-rules determination of 
quark masses has been incremental as has been 
the case for LQCD. Over time the perturba- 
tion expansion for the 2-point hadronic correc- 
tion functions has been carried out to higher or- 

(3) 

der, along with a better determination of &q Cd . 
Models for the hadronic spectral functions have 
been improved. The main limitation continues 
to be the lack of experimental data, with the 
one exception of r decays. Thus, one has to 
model the spectral function rather than measure 
it. This slow but steady progress was, as summa- 
rized in Table |[ reaching a consensus by 1996 at 
7n(2GeV) » 5MeV and TO s (2GeV) « 140MeV. 

Sum rule calculations proceed in one of two 
ways, (i) Using axial or vector current Ward 
identities one writes a relation between two 2- 
point correlation functions, where the constant 
of proportionality are the quark masses H|| . (ii) 
Evaluate a given correlation function both by 
saturating with known hadronic states and by 
evaluating it in perturbative QCD (PQCD) |g5[ . 
The PQCD expression depends on quark masses, 
and defines the scheme in which they are mea- 
sured. Systematic errors arise from the (i) finite 
order calculation of PQCD expressions, (ii) the 
scale fi above which perturbative and hadronic 
solutions are valid and can be matched on the 
average (duality), and (iii) the ansatz for the 
hadronic spectral function. 



Table 2 

Values and bounds on to and to, 



from sumrulcs. 



reference 



to (MeV) to s (MeV) 



1989 
•1 1995 
1995 
1995 
1996 
1997 
1997 
1997 
4| 1997 
42] 1997 
42] 1997 



= 6.2(0.4) 
= 4.7(1.0) 
= 5.1(0.7) 



= 4.9(1.9) 
>3.8-6 

> 3.4 
> 4.1-4.4 



= 138(8) 

= 144(21) 
= 137(23) 
= 148(15) 
= 91 - 116 
= 115(22) 

> 118-189 
> 88(9) 

> 104- 116 



In the last year, with the calculation of 
terms in the perturbative expansions, the value 
of Agpp settling around 380MeV, and a critical 
reappraisal of the systematic errors in the sum- 
rule calculations |3J], there has been a flurry of 
activity as shown in Table |[ The highlights of 
the new works are as follows. 

The calculation of 0(a^) terms and a detailed 
analysis of the convergence of the perturbation 
expansion suggests that the associated error is 
under control at « 10% level for \i > 2 GeV @. 

Colangelo et al. ]38| have extended the analy- 
sis of m s in p|j37]] by constructing the hadronic 
spectral function from known phase shift data. 
Similarly, Jamin ]39| ] has also used a different 
parametrization of the hadronic spectral function 
using this phase shift data. In both cases the re- 
analysis lowers the estimate of the strange quark 
mass significantly as shown in Table |2|. 

Prades ]^0) has repeated the analysis of to in- 
corporating the a^. corrections and using J^q]j D = 
380 MeV. He reports a slightly higher value than 
in H . This is because Prades chooses the duality 
point at /i 2 = 2GeV 2 , where to has a maximum. 
There is a significant decrease with /z, the number 
dropping to 4.3(1.7) at /z 2 = 3GeV 2 , and 3.8(1.5) 
at p? = 4GeV 2 . The rationale for the low choice 
of n is that contributions not included in the spec- 
tral function will bolster the answer for larger [i. 
This assumption needs to be substantiated. 

Finally, a number of calculations have used the 
positivity of the spectral function to derive lower 
bounds |3^,p4|-|43[] which depend on [i. Of these, 
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the most stringent were reported by Lellouch at 
this conference [^3|,0. These bounds rule out 
the 1996 n f = 2 lattice results for fj, < 2.8 GeV. 
The hard to resolve question is - what is the scale 
fi at which PQCD, and thus the bound, becomes 
reliable? Unfortunately, this question cannot be 
answered at present. 

8. Bottom quark mass 

There are two determinations of mj,. The 
NRQCD collaboration determine it from the 
Upsilon binding energy and the APE collabora- 
tion [|6| use HQET. The two results agree: 

mp(mp) = 4. 15(5) (20) GeV (APE) 

= 4.16(15) GeV (NRQCD) (27) 

The nice features of the NRQCD method are 

(i) determination of the scale from the spin av- 
eraged IP — IS splitting. Data show that these 
splittings are independent of the precise tuning of 
the input mass mP Hl and of the light quark action. 

(ii) tuning of m° H using the kinetic mass 

M =Yan{d 2 E/dp 2 y 1 . (28) 
The heavy quark mass is then defined in two ways 
m p B ole = 7^[m HlI - (E sim - 2E )] 

m MS = Z pole^MS m Pole 

mW = Z m m° H . (29) 

The drawbacks are that ZP ole ~* m , Z m , and E 
(the energy of a zero-momentum quark state) are 
only known in PQCD to 1-loop. The two meth- 
ods give consistent results, however, it is essential 
that the calculation be done at other values of (3, 
to test for stability under variations of a. 

The APE collaboration uses HQET to define 

I \ ( u pa a s(a)X\ ( a s (m b )\ 
m b (m b ) = I Mb—c H J II — I .(30) 

They explain how the linear divergence in the 
energy measured on the lattice, £, is cancelled 
by that in the perturbative series whose first 
term is a s (a)X/a. The residual renormalon am- 
biguity is cancelled by a similar one in the sec- 
ond factor that relates the pole mass to the MS 



mass. However, to take into account effects of 
higher orders in perturbation theory, they assign 
a systematic error of 200MeV to the final re- 
sult. We believe that their analysis of the can- 
cellation of the renormalon ambiguity also ap- 
plies to the NRQCD analysis, and with a similar 
residual uncertainty. The last issue concerning 
APE data is the stability under variations of the 
lattice spacing. We feel that by averaging over 
data at j3 = 6.0 — 6.4, which shows a marginally 
significant variation, the APE collaboration may 
have missed an equally important source of un- 
certainty. 

9. Charm quark mass 

There are two new determinations of the charm 
quark mass by the APE jlTj and the Fermilab 
collaborations E^j. The extraction by APE uses 
the same method as in |l6| . The new features are 
that they use both the WI and the HS methods 
and use non-perturbative estimates of Z's in the 
WI method. The drawback once again is that 
they average the data at (5 = 6.0, 6.2, and are not 
able to resolve discretization errors. 

The Fermilab collaboration uses a combination 
of the NRQCD and Fermilab approaches. In their 
TI clover results, they include 0(ma) corrections 
in the determination of 1-loop Z's, and calculate 
the q* in the matching factor using the Brodsky- 
Lepage-Mackenzie prescription, which they argue 
gets rid off possible infrared scales seen in the 
connection to pole mass. They find agreement 
in the extrapolated value between the two ways 
of calculating m c , and also argue that quenching 
errors are expected to be small, unlike in the case 
of light quarks. 

The final results are 

mf^imf 3 ) = 1.525(40) (125)GeV APE . 

m|p(m^) = 1.33(8)GeV Fermilab . (31) 

We do not consider the difference significant as it 
could easily be due to the different ways of setting 
the scale and/or due to the discretization errors 
that the APE data does not resolve. 

Finally, we would like to mention the work of 
Bochkarev and Forcrand j4Sj . They calculate m c 
by evaluating the correlation functions arising in 
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QCD sum-rules on the lattice. Their estimate, in 
the quenched approximation, is 



1.22(5)GeV 



(32) 



The errors include effects of discretization and fi- 
nite volume, and in fact are dominated by the 
scale uncertainty. Given that these results were 
obtained on relatively small lattices (16 3 x 32) and 
with low statistics (20 configurations), the quoted 
accuracy is impressive. This approach should be 
investigated further. 

10. Conclusions and Acknowledgements 

Reliable estimates of quark masses have two 
immediate phcnomcnological consequences, (i) 
Standard Model predictions of e'/e are very sen- 
sitive to m s + md J49|, and (ii) they constrain 
Supersymmetric models |5(J . ft is therefore excit- 
ing to report that the range listed in the Particle 
Data Book |5l| has been significantly reduced. 

In the last year the quenched estimates have 
been significantly improved, mainly due to the 
factor of 100 better data by the CP-PACS collab- 
oration. There has been progress by the APE and 
ALPHA collaborations in the non-perturbative 
determination of renormalization constants. The 
SESAM analysis has lead to a deeper understand- 
ing of how to extract quark masses from Nf = 2 
simulations, however, we are no further along in 
obtaining continuum limit estimates. Since the 
quenching errors are the least well understood, it 
is this issue that needs maximum attention. 

There has been progress in understanding er- 
rors in heavy quark analysis by the APE, Fermi- 
lab, and NRQCD collaborations. We now have 
"first generation" estimates of m c and mt>. 

An area that has seen no progress is the 
study of isospin breaking effects. Knowing 
whether m u = is important since a zero value 
solves the strong CP problem. The exploratory 
quenched calculations by Duncan et al p2[ , 
where an electromagnetic field was added to 
quantify isospin breaking, reported a non-zero 
value of m u : (m,d — m u )/m s = 0.0249(3) and 
m u /m,i = 0.512(6). Unfortunately, quenched 
calculations cannot address this question. The 
subtle point is that quark masses calculated from 



low energy phenomenology (or LQCD calcula- 
tions with chiral extrapolations based on trun- 
cated xPT expansions) include instanton induced 
finite renormalizations ]52[| . What one wants is 
m u defined at a high scale. LQCD can directly 
probe this, but only if simulations are done in- 
cluding very light dynamical fermions whereby 
the influence of instanton zero modes on the light 
quark propagation is properly included. 

Finally, it is exciting to see the rivalry between 
the LQCD and sum-rules estimates of quark 
masses heating up. We feel that both sides have 
made big strides in understanding and address- 
ing the various sources of systematic errors, and 
estimates of masses have been tightened. 

We thank C. Allton, L. Giusti, S. Gottlieb, 
S. Hashimoto, H. Hoeber, R. Kenway, A. Kron- 
feld, M. Luscher, S. Ohta, K. Schilling, G. Schier- 
holz, and T. Yoshie, for providing details of their 
data in advance. We also acknowledge informa- 
tive conversations with L. Lellouch, M. Luscher 
and S. Sharpe. We thank DOE and the ACL for 
support of our work. 
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